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Asymptotic Energy Decay in Inelastic Fluids 

J.A.G. Orza, R. Brito and M.H. Ernst 

Dpto. de Ffsica Aplicada I, Universidad Complutense 

28040 Madrid, Spain. 

The goal of the present publication is the comparison 
of the two existing theoretical predictions |[]||] for the 
long time behavior of the total energy E(t) in a large sys- 
tem of freely evolving inelastic hard spheres (IHS) with 
computer simulations of N = 50 000 hard disks. 

(i) The first theory, a mode coupling theory by Brito 
and Ernst 0], predicts that E ~ AT~ d / 2 for d > 2, where 
r(t) is the average number of collisions suffered by a par- 
ticle in time t (r is a nonlinear function of t, whose ana- 
lytic form is unknown at large t) and A is a known coef- 
ficient, which depends on the coefficient of restitution, r, 
and on the density, <j>. 
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FIG. 1. Energy vs r at packing fractions = 0.05 (up- 
per) and 0.4 (lower) at different inelasticities. The dotted 
lines represent the short time behavior (Haff's law), and the 
dashed-dotted lines the theory of [1]. 

Computer simulations of IiiE(t)/E(0) are shown in 
Fig. 1, where we have chosen units such that E(0) = 1. 
These plots confirm that the small fluctuation theory of 
jl| gives at fixed <j> quantitative predictions inside an 
r— window of width Ar ~ 0.1, centered around r ~ 
0.6, 0.75, 0.8, 0.85 at densities cj> ~ 0.05, 0.11, 0.245, 0.4 
respectively. At larger/smaller r— values, E(r) decays 
faster/slower than the prediction of |Q. The location tq 
is determined by two balancing effects: sufficiently high 
density to suppress large relative density fluctuations 
which increase the mean overall collision frequency dr / dt 
(compression of r— axis), competing with sufficiently high 
inelasticity, which favors local inelastic collapses (a finite 
number of collisions in an infinitesimal time). Such colli- 
sions merely increase the r— value (stretching of r— axis) 



at fixed t without advancing the N particle dynamics. 
At small density the relative density fluctuations become 
very large, as the clusters keep growing and compactify- 
ing, which invalidates the linear theory of ]lj. 

(ii) Ben-Nairn et al. show that the behavior of a 
fluid of inelastic hard rods is described by the totally 
inelastic (sticky) gas. Consequently, the energy decays 
at long times as i~ 2 / 3 in 1-D, and is independent of the 
coefficient of restitution. Moreover, they conjecture that 
the total energy E for 2 < d < 4 decays as Bt~ d / 2 with an 
unknown coefficient B, independent of inelasticity. The 
results of Q are only valid at asymptotically large times. 

Our goal is to test this conjecture against simulations 
of inelastic hard disks. The energy decay when plotted 
as log-E versus logt (see Fig. 2), gives the misleading im- 
pression that our simulations have reached their asymp- 
totic time dependence, and suggests that E ~ Ct~ a de- 
cays algebraically with a density-dependent exponent a, 
but with a and C independent of the dissipation, possi- 
bly corresponding to a sticky hard sphere fluid. At small 
densities (Fig. 2a), a ~ 1, which offers partial support 
for the conjecture, while at higher densities (Fig. 2b) the 
analysis seems to show a smaller exponent. 
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FIG. 2. Energy vs t at packing fractions (j> — 0.05 (up- 
per) and 0.4 (lower) at different restitution coefficient r. The 
dotted lines represent the short time behavior (Haff's law). 
Long-dashed line with slope —1 is plotted for reference. 

However, a more sensitive test is to plot E(0)/E(t) 
versus t, and test whether the curves for different values 
of the restitution coefficient r become linear in t, and tend 
to coincide for large times, i.e. become independent of r. 
This is done in Fig. 3 for three different packing fractions 
4> = 0.05,0.245,0.4, each for a range of r-values over a 
long time interval far beyond where Haff's law is valid. 
The results at low density can hardly be considered as 
evidence for the conjecture. The curves at intermediate 
density, = 0.245 show behavior that might be linear 



1 



in t, but depends strongly on the degree of inelasticity. 
Simulations at density <j> = 0.11 show behavior similar 
to the ones for <fi — 0.245. The behavior at the highest 
density at (f> = 0.4, in the time interval considered, looks 
roughly independent of the degree of inelasticity, but the 
curves show a tendency to diverge at later times. 
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FIG. 3. Long time behavior of the inverse of the energy vs 
t at packing fractions (j> = 0.05 (upper), <f) = 0.245 (middle) 
and <f) = 0.4 (lower) at different dissipations r. Short time 
behavior, 1/E(t) ~ (1 + t/t e ) 2 , as described by Haff's law is 
still visible. 



asymptotically large times in small systems. In thermo- 
dynamically large systems this crossover is never reached. 

In small systems we have observed that the energy de- 
cays as, E(t) ~ C(N,r)/t 2 for t — > oo, with a coefficient 
C that depends on N and r, and that differs in general 
from Haff's law. The law E(t) ~ t~ 2 seems to hold for 
small systems in any dimension larger than 1 (see also 
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In general, the curves in Fig. 3 do not seem to have 
reached their asymptotic form, and are not conclusive 
enough to support or refute conjectures about the asymp- 
totic ^-dependence being independent of the degree of 
inelasticity. 

Finally we observe that at asymptotically large times 
there is an important distinction between small and ther- 
modynamically large systems, for the following reason. 
The growth of patterns, i.e. vortices and density clus- 
ters, is controlled by diffusive modes [0, and typical di- 
ameters, L v (t) of these patterns grow as y/r. As soon 
as the system size L ~ L v (t), patterns start to interfere 
with their periodic images, and there occurs crossover to 
a steady state, which is fully determined by the (unphys- 
ical) periodic bounday conditions. This is the case for 
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